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Abstract. Results of simultaneous INTEGRAL and optical observations of the galactic microquasar SS433 in May 
2003 and INTEGRAL /RXTE observations in March 2004 are presented. Persistent precessional variability with a 
maximum to minimum uneclipsed hard X-ray flux ratio of ~ 4 is discovered. The 18-60 keV X-ray eclipse is found 
to be in phase with optical and near infrared eclipses. The orbital eclipse observed by INTEGRAL in May 2003 
is at least two times deeper and apparently wider than in the soft X-ray band. The broadband 2-100 keV X-ray 
spectrum simultaneously detected by RXTE/ INTEGRAL in March 2004 can be explained by bremsstrahlung 
emission from optically thin thermal plasma with kT ~ 30 keV. Optical spectroscopy with the 6-m SAO BTA 
telescope confirmed the optical companion to be an A5-A7 supergiant. For the first time, spectorscopic indications 
of a strong heating effect in the optical star atmosphere are found. The measurements of absorption lines which 
are presumably formed on the non-illuminated side of the supergiant yield its radial velocity semi-amplitude 
K v = 132 ± 9 km/s. The analysis of the observed hard X-ray light curve and the eclipse duration, combined with 
the spectroscopically determined optical star radial velocity corrected for the strong heating effect, allows us to 
model SS433 as a massive X-ray binary. Assuming that the hard X-ray source in SS433 is eclipsed by the donor 
star that exactly fills its Roche lobe, the masses of the optical and compact components in SS433 are suggested 
to be M v ~ 30Mq and M x w 9Mq, respectively. This provides further evidence that SS433 is a massive binary 
system with supercritical accretion onto a black hole. 
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1. Introduction 

SS433 is a massive eclipsing X-ray binary system at an 
advanced evolutionary stage. It is recognized as a super- 
critically accreting microquasar with a precessing accre- 
tion disk and mildly relativistic (v w 0.26 c) jets. Since 
its discovery in 1978 (Clark and Murdin 1978, Margon 
et al. 1979), this unique X-ray binary has been deeply 
investigated in optical, radio, and X-rays (see reviews by 
Margon 1984, Cherepashchuk 1988, 2002 and Fabrika 2004 
for more detail and references). 
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Among a dozen known microquasars, SS433 is distin- 
guished by its unique properties. 

1) In this system the optical star fills its Roche lobe 
and mass transfer occurs on thermal time scale with the 
huge rate of M ~ 10 _4 M Q per year and accretion onto the 
relativistic object is supercritical (Shakura and Sunyaev 
1973). 

2) Two strongly collimated (~ 1°) oppositely directed 
jets emanate with a velocity of 80,000 ± 1,000 km per 
second from the center of the accretion disk. The jets are 
observed not sporadically, as in most microquasars, but 
persist over tens of years with virtually constant velocity. 
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3) The accretion disk and relativistic jets regularly pre- 
cess with a period of 162.5 days. Both the precession pe- 
riod and the disk inclination angle to the orbital plane 
(~ 20°) on average remain constant over tens of years. 

4) The source exhibits orbital eclipsing periodicity 
with period P orb — 13 d . 082. The shape of the optical 
light curve varies significantly with precessional phase 
(Goranskij et al. 1998ab, Cherepashchuk and Yarikov 
1991). The orbital period is found to be very stable over 
more than 25 years, despite the high mass exchange rate 
between the binary components and strong wind mass loss 
(v fa 2000 km/s, M ~ 10 _4 A/ Q /yr) from the supercritical 
accretion disk. 

The unsolved puzzles of SS433 still remaining to be 
solved include: (1) the nature of the relativistic object, (2) 
the mechanism of collimation and acceleration of matter 
in jets to the relativistic velocity ~ 80, 000 km/s, (3) the 
nature of the precessional phenomena in this X-ray binary 
system (Chakrabarti, 2002). 

The INTEGRAL observations of SS433 in 2003 dis- 
covered a hard (up to 100 keV) X-ray spectrum in this 
supercritically accreting microquasar (Cherepashchuk et 
al. 2003, 2004), suggesting the presence of an extended 
hot (with a temperature up to 10 s K) region in the cen- 
tral parts of the accretion disk. These new data made it 
possible to compare the eclipse characteristics of SS433 at 
different energies: soft X-rays (2-10 keV, the ASCA data), 
soft and medium X-ray (1-27 keV, the Ginga data), hard 
X-rays (20-70 keV, the INTEGRAL data), and the opti- 
cal. This comparison allows us to investigate the inner- 
most structure of the supercritical accretion disk and to 
constrain the basic parameters of the binary system . 

This paper is organized as follows: Section 2 lists 
the participants of the observational campaign. Section 
3 describes the X-ray spectra and light curves obtained. 
Section 4 presents the analysis of hard X-ray precessional 
and orbital variability of SS433. Section 5 reports the re- 
sults of optical and infrared photometry. Section 6 de- 
scribes in detail the simultaneous optical spectroscopy of 
SS433 obtained with the 6-m BTA telescope and the mea- 
surements of the optical radial velocity curve. The impli- 
cations of our observations for binary masses in SS433 are 
summarized in Section 7, and Sections 8 and 9 give the 
discussion and conclusions, respectively. 



2. Observational campaign 

The coordinated multiwavelength observational campaign 
of SS433 was organized during the INTEGRAL observa- 
tions of the SS433 field in March-June 2003 with the par- 
ticipation of the following research teams. 

1. Special Astrophysical Observatory of the Russian 
Academy of Science (SAO RAS). High signal-to- noise op- 
tical spectroscopy of SS433 at the 6-m telescope. 

2. Crimean Station of the Sternberg Astronomical 
Institute. BVR-photometry at the 0.6-m telescope. 
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Fig. 1. IBIS/ISGRI spectrum of SS433 collected over 
INTEGRAL orbits 67-69 in May 2003. The best fit (solid 
line) is for an exponential cut-off with E c — 14 ± 2 
keV. Reduced \ 2 — 0-32 for 7 degrees of freedom (dof), 
CL ss 95%. 



3. Kazan State University. Optical photometry at 
the Russian- Turkish 1.5-m telescope RTT-150 at the 
TUBITAK National Observatory (Turkey). 

4. Infrared K-photometry was performed with the 1.1- 
m AZT-24 telescope of the Pulkovo Observatory in Italy, 
Osservatorio di Campo Imperatore. 

5. Radio monitoring at cm wavelength: the RATAN- 
600 radio telescope of SAO RAS. The source was found to 
be in its non-flaring (quiet) state in May 2003 (Trushkin 
2003). 

SS433 was within the field of view of the IBIS telescope 
during the INTEGRAL observations (March 2004) of the 
Sagittarius arm tangent region centered on I = +40°, b = 
0°. We included part of these observations in our analysis. 
In contrast to the observations in May 2003, radio moni- 
toring showed that SS433 had been in its flaring (active) 
state (Trushkin 2004). 

3. IBIS/ISGRI light curves and spectra 

The INTEGRAL IBIS/ISGRI data were analyzed us- 
ing both publicly available ISDC software (OSA-3 ver- 
sion) and the original software elaborated by the IKI 
INTEGRAL team (see Revnivtsev et al. 2004b for more 
detail). The latter allows better sensitivity in constructing 
maps and spectra. The main results can be summarized 
as follows. 

A. The IBIS/ISGRI spectrum of SS433 (25-100 keV) 
in May 2003 as processed by the OSA-3 software can be 
fitted by a power-law ~ E~ a with photon index a « 2.7± 
0.13, taking into account the 10% systematic errors. A 
satisfactory fit is also obtained using an exponential cut- 
off esqp(-E/E c ) with E c ~ 12 - 16 keV (Fig.QJ. 

The integrated hard X-ray luminosity in May 2003 was 
L x (18 - 60keV) - 4 x 10 35 erg/s, L x (60 - 120keV) - 
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Fig. 2. RXTE PCA and IBIS/ISGRI spectrum of SS433 
obtained during simultaneous RXTE/INTEGRAL ob- 
servations of SS433 in March 24-27, 2004. The best fit 
(solid line) is for bremsstrahlung emission from optically 
thin plasma with kT = 29.0 ± Q.lkeV (model "bremss" 
from the XSPEC package). Reduced \ 2 = 0.77 for 50 dof, 
CL « 88%. 
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Fig. 3. IBIS/ISGRI 20-40 keV (upper panel) and 40-70 
keV (bottom panel) count rates of SS433 (without back- 
ground subtraction) obtained in May 2003. The time ori- 
gin is MJD= 52763.95. The egress part of the X-ray eclipse 
was kindly provided by Dr. Diana Hannikainen. 

Table 1. Best- fit parameters of the joint 
RXTE+INTEGRAL spectrum of SS433 observed in 
March 24-27, 2004. The model is bremsstrahlung emis- 
sion and a Gaussian emission line (from the XSPEC 
package). The energy of the emission line was a free pa- 
rameter. The physical nature of the line is not discussed 
in the present paper. 



Parameter 


kT 


Emission line 


Line 






energy 


width 




(kev) 


(keV) 


(keV) 


Value 


28.9 ±0.7 


7.26 ±0.03 


1.2 ± 0.04 



2 x 10 35 erg/s (assuming the 5 kpc distance to SS433), 
which is about 10% of the soft X-ray jet luminosity. 

Using the OSA-3 software, we could not significantly 
detect the source in the JEM-X data. Instead, we made 
use of the more detailed RXTE observations of SS433 per- 
formed simultaneously with INTEGRAL in March 2004 to 
obtain the broadband 2-100 keV spectrum of SS433. The 
source was observed a few days after the disk maximum 
opening phase. The OSA-3 software also did not allow us 
to significantly detect the source at energies above 70 keV, 
so we used the IKI software to obtain the SS433 spectrum 
up to 100 keV. This software has proven its quality and 
efficiency in processing INTEGRAL observations of the 
Galactic center (Revnivtsev et al. 2004b). 

The resulting broadband X-ray spectrum of SS433 in 
March 2004 is shown in Fig. It can be adequately fitted 
by the bremsstrahlung emission from optically thin ther- 
mal plasma with kT ~ 30 keV ( model "bremss" from 
XSPEC package was used; see Table 0). The reduced chi- 
square value for 50 dof is ~ 0.8, which corresponds to a 
null hypothesis probability of ~ 0.9. 

B. The IBIS/ISGRI count rates of SS433 in May 2003 
are presented in Fig. [3] The X-ray eclipse at hard ener- 
gies is observed to be slightly narrower than the optical 
one, slightly broader than in the 4.6-27 keV Ginga energy 
range and displays extended wings (Fig. EJ) . This is oppo- 
site to what is found in ordinary eclipsing X-ray binaries 
(like Cen X-3, Vela X-l etc.), in which the X-ray eclipse 
duration decreases with energy. 

C. The eclipse depth is observed to be at least 80% in 
hard X-rays compared to ~ 50% in the 4.6-27 keV band 
(Fig. EP. 

D. The 25-50 keV X-ray flux increases from ~ 5 to ~ 20 
mCrab during March-May 2003 when the source precessed 
from the crossover (T2) phase to the maximum opening 
disk phase (T3). This modulation is ~ 2 times larger than 
observed in the 2-10 keV energy band (see Fig- EJ> - Thus, 
both precessional and eclipsing hard X-ray variabilities 
in SS433 exceed by ~ 2 times those in the standard X- 
ray band (Fig. |5J). This suggests a more compact vertical 
structure of the hard X-ray emitting region in the central 
parts of the accretion disk. 

4. Analysis of precessional and orbital variability 
of SS433 in the 25-50 keV energy range 

4.1. Precession variability 

As mentioned above, the 25-50 keV flux of SS433 varies 
from F m i n ~ 5 mCrab to F max ~ 20 mCrab, with 
F m i n and F max corresponding to precession phases of 
the disk seen edge-on (moment T2) and at maximum 
open (moment T3), respectively, yielding the ratio A pr = 
F m ax /F m i n ~ 4. In Fig. 0] we present all our available 
INTEGRAL observations of SS433. The left panel shows 
the 2003 data (INTEGRAL orbits 49, 51, 53, 56, 58, 60, 
67-70). The right panel shows the 2004 data (March 24- 
25, orbits 176 and 177). The uneclipsed flux of SS433 in 
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Table 2. Change of the precession amplitude of SS433 
A pr with energy 
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Fig. 4. Precessional hard X-ray (25-50 keV) variability of 
SS433. Left panel: March-May 2003. Right panel: March 
2004. The filled squares mark the T3 precession phase 
according to ephemeris by Goranskij et al. [1998ab]. The 
upper axis shows the precessional phase according to the 
same ephemeris. 
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Fig. 5. The precession and eclipse amplitudes of SS433 at 
different energies in the quiet state: A pr (squares and solid 
line) and A ec \ (triangles and dashed line). 



2004 (when the source was observed near the T3 phase) 
is comparable within the errors with the maximum flux 
observed near the T3 phase in 2003 (these phases are in- 
dicated by the filled squares). This suggests that the pre- 
cessional hard X-ray variability of SS433 stays constant at 
least over several precessional cycles. 

The amplitude of the precession variability in SS433 
in different X-ray bands appears to be monotonically in- 
creasing with energy (Fig. the squares and the solid 
line; Table EJ). For example, archive RXTE/ASM data 
collected over several years (only quiet states of SS433 
have been selected) indicate that A w {\.h — 5keV) = 1.6, 
A pr (5 - 12keV) = 2.4 (Fabrika et al. 2004). The errors are 
due to the SS433 flux at the crossover being heavily ab- 
sorbed (especially at the first crossover). The Ginga data 
(Kawai et al. 1989, Yuan et al. 1995) also fit this depen- 
dence (Table EJ. The monotonic increase of the precession 
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[1] - Fabrika et al. 2004; [2] - Kawai et al. 1989; [3] - Yuan et 
al. 1995; [4] - Cherepashchuk et al. 2003 



amplitude with energy is consistent with the model of an 
emitting region in SS433 as a cooling outflow. 

For a jet of finite length or consisting of individ- 
ual evolving fragments with finite life times (Lind and 
Blandford 1985; Begelman et al. 1984; Panferov and 
Fabrika 1997), the intensity of radiation in the jet rest 
frame is J co = (1 + z) 2+a J b s , where J b s is the observed 
intensity of the jet emission, z — —8v/i> is negative for the 
approaching jet and a is the spectral index (I v oc i>~ a ). 
With appropriate values of the spectral index, the preces- 
sion amplitude due to the relativistic beaming only in the 
case where both jets are always seen and are not screened 
by the accretion disk would be ~ 1.1, and in the case 
that the receding jet is fully screened by the disk it would 
be ~ 1.5. So the large (~ 4 times) precession variability 
found in the hard X-ray band 25-50 keV cannot be ex- 
plained only by the jet relativistic beaming and should 
be almost fully caused by geometric screening of the in- 
ner disk structure producing hard X-rays by the edge of 
the precessing thick accretion disk. The precession light 
curve in the 25-50 keV energy band exhibits a significant 
irregular variability which is apparently related to the non- 
stationarity of the outer edges of the geometrically thick 
accretion disk. 

4.2. Hard X-ray eclipse 

The INTEGRAL observations in May 9-11, 2003, enable 
us to study in detail its hard X-ray eclipse. We have used 
our data of May 9 and 11. Data for May 10, 2003, were 
kindly provided by Dr. Diana Hannikainen. The main 
eclipse falls on the precession phase ip pr = 0.09 (according 
to the ephemeris by Goranskij et al. 1998a), fairly close to 
the T3 (disk maximum opening) phase. 

The primary X-ray eclipse in the 18-60 keV energy 
band (IBIS/ISGRI data) is plotted in Fig. El The points 
are averaged over 5 and 10 INTEGRAL scientific windows 
(10000 and 20000 s, respectively). For comparison, we plot 
the 4.6-27 keV X-ray eclipse observed by Ginga (Kawai et 
al. 1989, Yuan et al. 1995) at about the same precession 
phase. The upper panel in Fig. H3 also shows the optical 
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Table 3. Primary X-ray eclipse depth A ec i variations with 
energy 
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1.0-4.6 


2.8 


1.7 


Ginga 


[1,2] 
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Ginga 
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2.9 


Ginga 


[1,2] 
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INTEGRAL 


[3] 



[1] - Kawai et al. 1989; 
Cherepashchuk et al. 2003 



[2] - Yuan et al. 1995; [3] 



(V) light curve observed simultaneously with INTEGRAL 
in Crimea and SAO (see below) . The upper panel demon- 
strates the absence of any noticeable lag between optical 
and hard X-ray primary eclipses. 

Fig. El shows that the INTEGRAL eclipse is up to two 
times deeper than that observed by Ginga in softer X-ray 
bands. The ascending branch of the hard 18-60 keV eclipse 
clearly shows some irregularities. These features are sim- 
ilar to those observed by Ginga in the softer 4.6-27 keV 
X-ray band. If this is indeed the case, we can take the up- 
per envelope of the ascending branch of the X-ray eclipse 
light curve as a representative shape of the eclipse by the 
optical star body. Note that the X-ray eclipse minimum 
observed by Ginga was consistent with the synchronously 
observed optical one (Aslanov et al. 1993). This proves 
that there is no appreciable phase shift between the mid- 
dle of the INTEGRAL and Ginga X-ray eclipses shown in 
Fig. El 

The X-ray eclipse depth as a function of energy is 
shown in Fig.[3]by the dotted line. The eclipse depth A ec i 
was defined as the ratio of the mean uneclipsed flux before 
the eclipse to the middle eclipse flux at the T3 precession 
phase of SS433 in the quiet state. These data are col- 
lected in Table El Note the tendency for A ec i to increase 
with energy, which may mean that the hotter parts of the 
X-ray emitting region located closer to its base are totally 
eclipsed by the optical star. 

The dependence of the precession amplitude and X- 
ray eclipse depth on energy is the same, suggesting similar 
screening efficiency of the X-ray emitting region by the ac- 
cretion disk edge and by the optical star. This means that 
the size of the accretion disk can be comparable with that 
of the optical star. Such a conclusion has been also in- 
ferred from optical photometry (Fabrika & Irsmambetova 
2002; Fabrika 2004). 

4.3. The geometric model of the system 

To interpret the X-ray light curve of SS433 we used a ge- 
ometric model applied earlier to the analysis of the Ginga 
X-ray eclipse (Antokhina et al. 1992). We consider a close 
binary system consisting of an opaque "normal" star lim- 
ited by the Roche equipotential surface and a relativistic 
object surrounded by an optically and geometrically thick 
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Fig. 6. The primary X-ray eclipse of SS433 in the I8- 
60 keV energy band. IBIS/ISGRI data (May 2003). 
Upper panel: X-ray light curve averaged over 20000 s 
(10 INTEGRAL science windows, SCW) superimposed on 
the simultaneously obtained V optical photometric light 
curve (Crimea, SAO). Bottom panel: the same hard X-ray 
eclipse light curve averaged over 10000 s (5 SCW) super- 
imposed on the Ginga 4.6-27 keV eclipse (filled circles; 
from Kawai et al. 1989, Yuan et al. 1995) taken at about 
the same precession phase. The INTEGRAL data is the 
same in both panels, but the averaging is different. 




Fig. 7. Geometrical model of the accretion disk and its 
"corona" . 

"accretion disk". The "accretion disk" includes the disk 
itself and an extended photosphere formed by the outflow- 
ing wind. The orbit is circular, the axial rotation of the 
normal star is synchronized with the orbital revolution. 

The disk is inclined with respect to the orbital plane by 
an angle 9. The opaque disk body (see Fig.0) is described 
by the radius and the opening angle lu. The central ob- 
ject is surrounded by a transparent homogeneously emit- 
ting spheroid with a visible radius Tj and height bj , which 
could be interpreted as a "corona" or a "jet" (without any 
relativistic motion). 

Here rj , bj and rd are dimensionless values expressed in 
units of the binary separation a. The radius of the normal 
(donor) star is determined by the relative Roche lobe size, 
i.e. by the mass ratio q = M x /M v (M x here denotes the 
mass of the relativistic object). 

Only the "corona" is assumed to emit in the hard X- 
ray band, while the star and disk eclipse it during orbital 
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and precessional motion. During precession the inclination 
of the disk with respect to the observer changes, causing 
different conditions of "corona" visibility. Observations of 
the precessional variability can thus be used to obtain a 
"vertical" scan of the emitting structure, restricting the 
parameters bj and u>. The orbital (eclipse) variability ob- 
servations scan the the emitting structure "horizontally" , 
restricting possible values of r<j, u>, q and rj. The joint 
analysis of precessional and eclipse variability enables us 
to reconstruct the spatial structure of the region in the 
accretion disk center where the hard X-rays are produced. 

The position of the components of the system relative 
to the observer is determined by the binary orbit incli- 
nation angle i = 78.8°, the disk inclination angle to the 
orbital plane 8 = 20.3°, the precession phase ip pr with 
ip pr = at the maximum disk opening of SS433 (T3, max- 
imum separation between the moving emission lines) and 
Vv = 0-34, 0.66 when the disk is seen edge-on (at the 
moving emission line crossover moments Tl and T2, re- 
spectively). 

To analyze the X-ray eclipse and precession light curve 
observed by INTEGRAL we have used the following three 
model shapes of the hard X-ray emitting region which can 
be parametrized by two parameters rj and bj : long narrow 
jets(rj <C bj ~ a), short narrow jets (rj <C bj <C a), and 
short thick "jets", or corona (rj > bj). For all models the 
formal \ 2 minimum is reached for a maximum allowed ac- 
cretion disk radius r ( i and the disk opening angle ui ~ 90° 
(see discussion below). The accretion disk radius was lim- 
ited by the distance from the relativistic object center to 
the inner Lagrangian point, which exceeds the radius of 
the Roche lobe around the compact star. This is natu- 
ral for a supercritical accretion disk with a strong radial 
outflow (Zwitter, Calvani & D'Odorico, 1991). Further ex- 
tension of the accretion disk is not justified because of the 
presence of narrow absorption lines of the optical A5-A7 
component in the spectrum, which means that there is no 
appreciable screening of the optical star by the disk. Such 
an extension of the disk would mean that the binary sys- 
tem is on the common envelope evolutionary stage, and so 
the narrow absorption lines from the optical star surface 
could not be observed. The orbital period stability over 30 
years of observations also suggests against the presence of 
a common envelope in SS 433. 

The results of our analysis can be summarized as fol- 
lows. 

1. Long narrow jets (rj <C bj ~ a) outflowing from 
the central object can be only partially screened by the 
disk during its eclipse and precession. Neither the observed 
eclipse depth nor the precession amplitude A pr ~ 4 can 
be explained by this model. (See Fig. |HJ . 

2. Short narrow jets (rj <C bj <C a) where bj is found 
from the condition that the precession variability ampli- 
tude is A pr — 4. The theoretical eclipse light curves for 
this model for different mass ratios q are displayed in Fig. 
El For any q this model is inconsistent with the observed 
18-60 keV X-ray eclipse shape. 




q=0.05 



Fig. 8. X-ray eclipse fit with a model of narrow long X- 
ray emitting jets. Precession variability cannot be repro- 
duced by this model although the eclipse shape is fitted. 
The point at the center of the eclipse was taken from the 
Galactic Plane Survey scan in the beginning of the 70th 
INTEGRAL orbit 




Fig. 9. The same as in FigO but for a model of nar- 
row short X-ray emitting jets for different values of q. 
Precession variability amplitude A pr = 4 is reproduced, 
but the eclipse shape is not. 



3. An extended corona (rj > bj) gives the best fit 
for the observed hard X-ray precession amplitude (A pr ~ 
4 ) and the observed eclipse shape (Fig. I10|) . Minimum 
residuals in this model are attained for the mass ratio q = 
m x /m v w 0.2 fFig.ITT]!. which is close to that derived from 
fitting the Ginga and ASCA observations with a narrow 
jet model (Kawai et al. 1989, Kotani et al. 1996). If we take 
into account only the upper envelope of the ascending X- 
ray eclipse branch (see above), the fit with q « 0.3 seems 
to be acceptable as well (Fig. an d 11 lf> . 

The parameters bj and to appear to be related and 
cannot be determined separately. Formal \ 2 decreases as 
u — > 90°, without any essential change of q for u ~ 80°. 
Obviously, the solution with u> = 90° (planar disk with 
zero thickness) could not provide any precessional vari- 
ability, so we accepted ui ~ 80°. This also agrees with the 
estimates of the thickness of accretion disk from the rapid 
optical and X-ray variability corellation (Revnivtsev et al., 
2004a). 

From this analysis we conclude that in the framework 
of our geometrical model of the hard X-ray emitting re- 
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Fig. 12. V-light curve of SS 433 obtained at the RTT150 
telescope (TUBITAK National Observatory, Turkey) si- 
multaneously with INTEGRAL observations. Bottom: 
photometry of control stars (Vn3 = 12.93, Vjv4 = 12.70). 



Fig. 10. The same as in Fig. |H1 but with a model of an 
extended oblate X-ray emitting region ( "corona" ) for dif- 
ferent q. The precession variability amplitude is A pr = 4. 
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Fig. 13. Rapid photometric variability of SS433 during 
the night 05/06 May 2003 



Fig. 11. x 2 map (26 dof) in the plane (rj,q) for models 
with an extended oblate X-ray emitting region. The outer 
edge of the accretion disk limited by the distance to the 
inner Lagrangian point is shown by the solid line. The 
precession variability amplitude is A pr = 4. 

gion in SS433, the observed precession amplitude and X- 
ray eclipse shape arc best reproduced by a broad oblate 
corona above an optically thick accretion disk. The best 
fit parameters are rj ~ 0.3, bj ~ 0.1 and uj ~ 80°, the 
mass ratio q — 0.2 — 0.3. 

Physically, such a corona could be thought of as hot 
rarefied plasma filling the funnel around the jets. In this 
picture outer parts of long thin jets that are not screened 
by the funnel walls are responsible for emission in soft X- 
rays, while the more extended hot corona inside the funnel 



produces hard X-ray flux. The estimation of the funnel size 
~ 10 12 cm was obtained by Revnivtsev et al. (2004a) from 
the analysis of rapid non-coherent variability in SS433. 
The opening angle of the funnel can be quite large, so 
the apparent size of the outer funnel that we model here 
as an extended oblate corona over an opaque disk can be 
comparable with the accretion disk size, as our best fit 
parameters suggest. The nature of such a plasma inside 
the funnel and how it maintains its stationarity requires 
further studies. We also note that actual detection accu- 
racy of the hard part of the broadband X-ray spectrum 
of SS433 obtained by INTEGRAL (Fig.H) is insufficient 
to definitely rule out the presence of the two-component 
structure of the X-ray emitting region inferred from our 
analysis. 
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Fig. 14. JHK-photometrical light curve of SS433 ob- 
tained by AZT-24 1.1-m IR telescope in July- August 
2003 (Campo Imperatore, Raly). The dashed vertical 
lines mark the primary eclipse of SS433 according to the 
ephemeris by Goranskii et al. (1998a). 

5. Optical and IR photometry 

5.1. Optical photometry 

Photometry of SS433 was performed simultaneously with 
INTEGRAL observations by the Russian-Turkish RTT- 
150 telescope of Kazan University at the TUBITAK 
National Observatory, Turkey. Observations were made 
using commercial camera (model DW436, www.andor- 
tech.com) which is termoelectrically cooled to —60° C. 
The CCD was provided for RTT150 by Max-Planck 
Institute for Astrophysics (Garching, Germany). It is a 
low noise back-illuminated model from EEV with 2048 x 
2048 pixels of 13.5 /j, size. The full field of view is 8 x 8 
arcmin with a frame reading time of 40 sec at 2 x 2 bin- 
ning. To increase the time resolution only parts of the field 
of view with reference stars N 1,2,3,4,5 ( Leibowits and 
Mendelson, 1982) around SS433 were captured. The ob- 
tained V-light curve is presented in Fig. El Strong (<~ 0.15 
mag) intranight variability of the source on timescales 
~ 100 s - 100 min is clearly detected (see Fig. ^|for the 
night 05/06 May 2003). 

The mean V light curve during the eclipse obtained 
at the Crimean Laboratory of SAI (Nauchnyi, Crimea) si- 
multaneously with the INTEGRAL observations is shown 
in the upper panel of Fig. [5J The detector was a pulse- 



counting, single-channel broad-band WBVR photoelec- 
tric photometer installed on a Zeiss-600 reflector. The 
photometer and software were designed and manufac- 
tured by I.I.Antokhin and V.G.Kornilov (SAI). The data 
were collected in the standard way using differential tech- 
niques. The main comparison star was CI with magnitude 
V = 11.51 (Gladyshev et al. 1980). The typical rms mea- 
surement errors in V were 0.02-0.03. The mean V light 
curve during the INTEGRAL observations corresponds to 
a quiet state of SS433. The optical eclipse minimum is ob- 
served at JD = 2452770.863, as predicted by the orbital 
ephemeris given by Goranskij et al. (1998ab). 

5.2. IR-photometry 

Near IR observations of SS 433 were obtained at AZT-24 
1.1m telescope in Campo Imperatore (Italy) during the 
period of July- August 2003. Some observational data were 
obtained also in October - November 2003. The AZT-24 
telescope at the Campo Imperatore Observatory located 
2150 m above sea level (a cooperation between Rome, 
Teramo and Pulkovo Observatories) is used for photo- 
metric studies of variable sources at near infrared (NIR) 
wavelengths. The telescope is equipped with a SWIRCAM 
NIR 256x256 pixels imaging camera mounted at the focal 
plane of AZT-24. This camera was built at the Infrared 
Laboratories in Tucson (Arizona, USA). The SWIRCAM 
is equipped with a PICNIC array, an upgrade of the 
NIGMOS detector, with a working range of 0.9 - 2.5 /i. 
It yelds a scale of 1.04 arcsec per pixel resulting in a 
field of view of 4x4 sq. arcmin. The observations were per- 
formed through standard Johnson JHK broadband filters. 
The NIR monitoring of SS433 started several orbits after 
the INTEGRAL observations. The JHK light curves are 
presented in Fig. ^] The observations were done close to 
the crossover precession phase, where the orbital modu- 
lation appears to be significantly reduced. Coincidence of 
the minima in the IR and optical light curves support our 
general geometrical model. 

6. Optical spectroscopy 

6.1. Spectral observations 

The optical spectroscopy of SS433 was performed at the 
6-m telescope of the Special Astrophysical Observatory 
during 6 nights on April 28 and May 9-13 2003 simulta- 
neously with the INTEGRAL observations at the preces- 
sional phase corresponding to the maximum disk open- 
ing angle (T3). The Long-Slit spectrograph (Afanasiev et 
al. 1995) at the telescope prime focus equipped with a 
1000x1000 Photometries CCD-detector was used to obtain 
spectra with a resolution of 3 A (1.2 A/pix). Standard 
techniques were used for spectral reduction and calibra- 
tion. 

Table 0] lists the spectral observations and includes 
date, JD of the middle of observation, orbital and proces- 
sional phases, number of spectra during the night and the 
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mean exposure, and spectral range used. We have taken 
spectra in the blue region as they are the most informa- 
tive when searching for the donor star absorption lines 
(Gies et al., 2002) and include the Hell A4686 emission 
line. On two nights, we obtained spectra in the red region 
to determine the orientation of relativistic jets and iden- 
tify the moving lines in all our spectra. The blue spectral 
range was shifted redward in the May 2003 observations 
because of the rising Moon. The signal-to- noise ratio in our 
spectra averaged over one night at A = 4250Ais » 60 per 
resolution element in May 10-13 and it is rs 80 on other 
dates. The signal-to-noise ratio increases towards longer 
wavelengths. 



6.2. Search for optical star lines 

The strongest lines in SS433, hydrogen, Hel, Hell and 
some Fell lines, appear in emission. Emission line widths 
in SS433 have FWHM ~ 1000 km/s, corresponding to 
~ 10 A in the blue region. The emission lines are formed 
in the accretion disk wind and in gas streams in the binary 
system. A detailed description of the SS433 spectrum can 
be found in Fabrika (2004). In precession phases close to 
the crossovers (T2), hydrogen (H/3 and upper lines), Hel 
and Fell lines have narrow blueshifted absorption compo- 
nents with P Cyg or shell-like profiles, which are formed in 
gas outflow from the accretion disk. Weaker Fell lines ap- 
pear as pure absorptions in these precession phases. The 
absorption components are also strengthened in orbital 
phase <j) ft 0.1; they are probably formed in the disk wind 
interacting with the star (Fabrika et al. 1997). 

The absorption components of the shell-like line pro- 
files and even the pure absorption lines are likely be 
formed in the disk wind too, which follows from the strong 
dependence of their intensities (Crampton & Hutchins, 
1981) and radial velocities (Fabrika et al. 1997) on the pre- 
cession phase. Therefore, a search for the donor star spec- 
trum in SS433 should be made with a certain caution. The 
search should be made among the weakest (photospheric) 
lines and in precession phases where the extended disk and 
its outflowing wind do not intersect the line of sight. In 
addition, stellar lines should be stronger in the middle of 
the accretion disk eclipse and weaker outside the eclipse 
(Gies et al. 2002). The main criterion that the absorption 
lines are formed in the stellar photosphere is a sign of the 
"correct" behavior of their intensities and radial velocities 
with orbital phase. 

Gies et al. (2002), Hillwig et al. (2004) detected weak 
absorption lines of Till, Fell, CrII, SrII, Cal, Fel in the 
blue region 4100-4600 A of the spectrum. They observed 
SS433 near the primary eclipse and in precession phases 
where the disk is the most " face-on" . These lines are be- 
lieved to belong to the donor star of SS433 with an A-type 
spectrum. 

One may also try to search for stellar lines in the yellow 
and red regions among the weakest metallic lines, however 
these regions are crowded with strong moving lines of hy- 



drogen and Hel formed in the relativistic jets. The moving 
lines are extremely variable and have structured profiles. 
This complicates the search for weak lines and the study of 
their behavior from date to date. For this reason we have 
investigated only the bluer parts of our spectra. The most 
informative region for stellar line studies is 4200-4340 A 
where no strong emission lines and no appreciable moving 
lines appear at these precession phases. 

The 4100-5300 A SS433 spectra have been com- 
pared with those of the known galactic supergiants 
with different temperatures. We used publicly available 



(http://webast.ast.obs-mip.fr/stelib/) library of stellar 
spectra (STELIB, Le Borgne et al. 2003). These spectra 
were taken with the same spectral resolution (3 A) as our 
spectra of SS433. For identification of lines we have used 
spectral atlases of Deneb (A2 Iae) (Albayrak et al. 2003) 
and o Pegasi (Al IV) (Gulliver et al. 2004). Spectra of 
SS433 summed over one night were used for the line iden- 
tification. 



7. Spectral type of the optical star 

Fig. E| shows spectra of four supergiants and two SS433 
spectra taken outside and inside the eclipse, each aver- 
aged over two nights. The eclipse-out spectrum was aver- 
aged over nights 09.05.2003 and 12.05.2003 corresponding 
to orbital phases 0.89 and 0.12 evenly spaced from the 
eclipse center (the top spectrum). The in-eclipse spectrum 
was averaged over nights 28.04.2003 and 11.05.2003 cor- 
responding to orbital phase 0.05 (spectrum second from 
the top in the figure). The supergiants shown in Fig. ^| 
are HD36673 (F0I, 7,200 K), HD39866 (A2I, 8,400 K), 
HD87737 (A0I, 9,700 K) and HD164353 (B5I, 22,000 K). 
All spectra were normalized and shifted along the vertical 
axis for better visualization. 

With our spectral resolution, practically all absorp- 
tion lines appear as complex blends. In Fig. E| the lines 
observed out of the eclipse are marked with short vertical 
bars; those seen in the eclipse or present in both spectra 
are marked with long bars. One can see that the strongest 
absorptions in the spectra of supergiants are also present 
in SS433 and these absorption lines are deeper in the spec- 
trum during the primary minimum. FeIIA4233 emission is 
seen (FWHM- 10 A), CIIA4267 and NIIIA4196, 4200 + 
HeIIA4200 emissions are also marginally present. 

The intensity of absorption lines during the eclipse al- 
lows us to estimate the effective temperature of the donor 
star to be T < 9000 K. This limitation was obtained 
from the fact that the absorption lines in SS433 cannot 
be deeper than those in spectra of supergiants, because 
the accretion disk is not totally eclipsed in the primary 
minima. The relative intensities of the strongest metallic 
absorption lines indicate a temperature of T = 8000 ± 500 
K, implying that the optical spectral class of the compan- 
ion is A5-A7I. 
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Table 4. Journal of spectral observations 



Date 


JD (mid) 


Orbital 


Precession 


N 


Exposure 


Range (A) 




2450000+ 


phase 


phase 




(sec) 




28.04.2003 


2758.491 


0.054 


0.973 


8 


900 


4160-5400 


09.05.2003 


2769.466 


0.893 


0.041 


3 


600 


5170-6410 


09.05.2003 


2769.507 


0.896 


0.041 


4 


900 


3970-5220 


10.05.2003 


2770.458 


0.969 


0.047 


4 


600 


5170-6410 


10.05.2003 


2770.501 


0.972 


0.047 


4 


900 


3970-5220 


11.05.2003 


2771.490 


0.048 


0.053 


8 


900 


4180-5430 


12.05.2003 


2772.508 


0.126 


0.060 


4 


900 


4180-5430 


13.05.2003 


2773.492 


0.201 


0.066 


4 


900 


4180-5430 
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Uaue 1 ength. A 

Fig. 15. Spectra of SS433 and spectra of four supergiants of known temperatures (Le Borgne et al., 2003) for a 
comparison. The top spectrum was taken outside the eclipse and averaged over nights 09.05.2003 and 12.05.2003 
(evenly spaced from the eclipse center orbital phases 0.89 and 0.12). The second from the top spectrum was taken 
inside the eclipse and averaged over nights 28.04.2003 and 11.05.2003 (the orbital phase 0.05). 



7.1. Evidence for a heating effect in the illuminated 
hemisphere 

The analysis of different absorption lines of the optical star 
reveals a strong heating effect in the illuminated star's at- 
mosphere. During the disk eclipse egress low-excitation 
absorption lines strongly weaken. The stellar hemisphere 
illuminated by the bright accretion disk in SS433 prob- 
ably has a temperature of ~ 20, 000 K, as the presence 



of CIIA4267 absorption + emission line suggests (the top 
spectrum in Fig. ll5fl . This absorption line is the strongest 
one in this spectral region among supergiants with tem- 
peratures T> 15, 000 K. The evolution of the blend A4273 
(FeIA4271.7 + FeIIA4273.3 + CrIIA4275.5) is seen: in the 
eclipse center the low-excitation Fel line is stronger, while 
out of the eclipse Fell+Crll lines are enhanced. Other Fel 
lines also appear only in the eclipse. 
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Fig. 16. Evolution of absorption lines in the optical spectra of SS433 taken at different orbital phases (see Table 
For comparison, spectra of two supergiants (8,400 K and 22,000 K) are shown in the bottom. Radial velocity of the 
optical star is traced by different absorption lines (vertical bars). On the left panel, the third from bottom is the SS433 
spectrum averaged over all nights with taking into account the radial velocity shifts from Fig. ^] See text for more 
details. 



These effects are illustrated by Fig. ^jj in which 
spectra of the standard stars with effective tempera- 
tures 8,400 K and 22,000 K are shown together with 
spectra of SS433 obtained on 9.05, 10.05, 28.04+11.05, 
12.05 and 13.05 (orbital eclipses fell on nights 10.05 and 
28.04+11.05). The orbital phase rises from bottom to top 
in the figure. There are emission lines in these spectral 
fragments — the broad structured emission of CIIA4267, 
and stronger emissions FeIIA5169+MgIA5167, 5173 and 
FeIIA5197. Note that emission lines, which presumably 
are formed in the disk wind, should move blueward (in 
phase with the compact object) in the orbital phases of 
our observations. In contrast, absorption lines, which are 
formed in the donor photosphere, must move red ward dur- 
ing our observations. The strength of the CIIA4267 ab- 
sorption component (marked by vertical bars in Fig. Iltif) 
increases notably out of the eclipse. The complex blend 
A4273 (Fel+Fell+Crll) is marked by the dotted lines in 
the figure. 

Fig. ^3 illustrates absorption lines (marked 
by the vertical bars) the blend A4314 of 

TiIIA4313+FeIIA4314+TiIIA4315; ScIIA4320 line and 
the blend ScIIA4325 + FeIA4326, 4327. However the last 



blend may be distorted by the strong blue wing of H7 
in the middle of eclipse. Other absorption lines shown in 
Fig. H3 (right panel) are A5154 (Till+Fell) and A5227 
(Till+Fel). It is seen that these lines shift redward with 
time. 

In Fig.^Dwe also show the spectrum of SS433 averaged 
over all nights of observations (the third from bottom). 
When averaging spectra obtained on individual nights we 
shifted them to zero velocity according to the radial veloc- 
ity curve from Fig. 1181 The average spectrum has a better 
signal-to-noise ratio. All the absorpton lines remain in the 
spectrum and they are similar to the corresponding lines 
in the spectrum of the standard supergiant (8,400 K). The 
shallow emission line (presumably CIIA4267) has an ab- 
sorption component. In our recent spectral observations 
of SS433, which were carried out around primary mimima 
on August 24, 2004 and September 6, 2004 (the preces- 
sion phase is the same, i.e. the maximum disk opening) 
we confirmed the presence of the absorption lines studied 
in this paper. Both the spectral resolution and signal-to- 
noise ratio are better in the new observations. Analysis of 
the new data will be published elsewhere. 
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Fig. 17. Radial velocity curve of the optical companion of 
SS433 obtained from 4 individual absorption lines over 6 
nights. 




Fig. 18. Mean radial velocity curve of the optical com- 
panion of SS433 measured from 22 individual absorption 
lines. The accretion disk radial velocities as measured 
by HcIIA4686 emission (Fabrika and Bychkova, 1990) is 
shown by the dotted line. 



7.2. Measurements of the radial velocity of the optical 
star 

It is important to stress that the strongest absorptions 
(Fell, Hel) are disfavored for the optical star radial veloc- 
ity analysis, as they most probably formed in the powerful 
outflowing disk wind (Fabrika et al. 1997). We studied ra- 
dial velocities of weak absorption lines of metals. It is easy 
to confuse the line identification when one traces a line 
from date to date because of the spectral variability and 
orbital motion of the companions. In Fig. 1171 we present 
individual absorption line radial velocities for four best- 
suited lines averaged over one night, as a function of the 
orbital phase for 6 nights. 

The radial velocity curve measured by the most reli- 
able collection of 22 absorption lines in the spectral range 
4200-5300 A is shown in Fig.UHl Note that with our spec- 



tral resolution the absorption lines in the SS433 spectrum 
are seen as blends containing 2-3 lines, and we measured 
their relative radial velocities because of unknown " labo- 
ratory" wavelengths of the blends. Relative intensities of 
the lines change from date to date through the eclipse. 
For this reason we measured only those lines and on only 
those dates where the lines were the most convincingly 
detected. The radial velocity curve (Fig. I18fl has been 
obtained by coadding radial velocity curves of individual 
lines, and the zero radial velocity of all lines on 11.05.2003 
(the orbital phase 0.048) was adopted. If a line was not 
detected on 11.05.2003, the zero radial velocity was as- 
sumed for 28.04.2003 as well (the orbital phase 0.054). If 
the measured radial velocity on 11.05.2003 disagreed with 
the whole radial velocity curve of a given line, an addi- 
tional shift for the curve was applied. However no shifts 
in excess of 20 km/s were done. This is the reason why the 
radial velocity error in 11.05.2003 is small, but not zero. 

Thus, the radial velocity curve for absorption lines in 
Fig. 1181 consists of individual radial velocity curves. Each 
point of the curve includes from 10 to 17 individual mea- 
surements. Most measurements were carried out at the 
primary minimum as there the absorption lines are deeper. 

The derived radial velocity semi-amplitude of the 
donor star is K v — 132 ± 9 km/s, the gamma-velocity 
of the binary system is w 7 = 14 km/s with a formal fit 
uncertainty of 2 km/s. The absorption line radial veloc- 
ity transition through the 7-velocity occurs at the mid- 
dle of the optical eclipse (0b — 0.07), confirming that the 
lines actually belong to the donor star. Note that in the 
gamma-velocity u 7 and in the transition phase 4>h some 
systematic errors can be present because of the method 
used. However they are less than 30 km/s (no bigger shifts 
were applied to radial velocity curves of individual lines) 
and 0.05 correspondingly. The main observational bias in 
studying radial velocities of the donor star is the faintness 
of absorption lines outside the primary minima and strong 
intrinsic spectral variability of SS433. To confirm the de- 
rived radial velocity curve further spectral observations 
are needed. 

Our results confirm the earlier determination of K v 
by Gies et al. (2002) which was carried out also at the 
maximum disk opening phases. Note that spectroscopic 
observations by Charles et al. (2004) were performed at 
the crossover phase of SS433 when the accretion disk is 
seen edge-on. Such a phase is disfavored for the donor star 
radial velocity analysis as strong gas outflows contaminate 
the disk plane; selective absorption in this gas affects the 
true radial velocity of the donor star. 

The heating effect of the donor star also distorts the 
radial velocity semi-amplitude. The analysis (Wade and 
Home 1988, Antokhina et al. 2005) indicates that the true 
value of the semi-amplitude of the radial velocity curve as 
derived from these absorption lines can be reduced to ~ 85 
km/s (see below). 

Recently Hillwig et al. (2004) obtained an estimate of 
the radial velocity amplitude of the donor star K v — 45 ± 
6 km/s and the gamma- velocity u 7 = 65 ± 3 km/s. In our 
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Fig. 19. The spectroscopic center of the normal star (cir- 
cle) heated by its companion is shifted with respect to the 
gravity center (cross). The arc arrows indicate the orbital 
motion. Zone "A" is the most effectively illuminated by 
the compact source, disk and jet. Zone "B" is also illumi- 
nated by outer parts of the disk and jet and by scattered 
radiation. The low excitation potential absorption lines 
used for radial velocity measurements can be formed in 
zone "C" which is not heated. Due to this effect, the spec- 
troscopically derived radial velocity of the optical star is 
overestimated. 

opinion, some strong lines (which actually are shell-like 
lines on the background of faint emission) could contribute 
to the radial velocity cross-correlations of Hillwig et al. 
(2004). For example, the two strongest lines in their Fig. 5 
FeIIA4550 and FeIIA4584 are in emission in our spectra in 
the middle of eclipse. 

The KPNO-2003 observations of Hillwig et al. (2004) 
were made in the orbital phases 0.85-1.07. Absorption lines 
formed in the rotating and extended envelope of the donor 
could be probably observed. In these orbital phases, the 
side of the envelope projected onto the strong continuum 
source (the accretion disk) moves away from the observer. 
Then the amplitude of the line shift can be expected to 
be ~ l/2V equa t ~ 60 km/s. This might explain the posi- 
tive shift of the system velocity +65 km/s found by those 
authors. 

8. Implications for binary masses 

Comparison of radial velocities of the accretion disk (K x — 
175 km/s, Fabrika and Bychkova 1990) and optical star 
(K v = 132 km/s) yields the mass ratio in the SS433 sys- 
tem q = m x /m v = K v /K x = 0.75 (here m x and m v stand 
for relativistic object and optical star mass, respectively). 
Taken at face value, these q, K v values would lead to the 
optical star mass function /„ « 3.12M© and the binary 
component masses m x — 18M Q , m v = 24M Q . However, 
such a large mass ratio is in a strong disagreement with 
the observed duration of the X-ray eclipse, which suggests 
a much smaller mass ratio q ~ 0.2 — 0.3. We stress that 
the binary inclination angle in SS433 (i = 78°. 8) is fixed 
from the analysis of moving emission lines. 

There are two possibilities: (1) either the model we 
used to fit X-ray eclipses should be modified, or (2) the 
value of K x and K v are influenced by additional phys- 
ical effects. Although there are some reasons to modify 



the model (e.g., the asymmetric shape of the hard X-ray 
eclipse, which may suggest an asymmetric wind outflow 
from the illuminated part of the optical star and wind- 
wind collision) , wc shall consider here only the second pos- 
sibility. The reason is that the actual value of K v should 
be decreased in order to account for the observed heating 
effect. 

Let us assume the mass ratio in the system to be 
q = 0.3, as with this value we can satisfactorily describe 
the width of X-ray eclipse and the observed X-ray preces- 
sion amplitude (see Fig.HTJI). Taking K v = 132 km/s yields 
/„ w 3.1M and m x = 62M , m v = 206M Q , K x = 440 
km/s. Clearly, this is an unacceptable model. Now let us 
decrease K v down to 85 km/s, the lower limit that follows 
from a more accurate treatment of the heating effect in 
the radial velocity curve analysis (Wade and Home 1988, 
Antokhina et al. 2005). This would yield a better fit with 
m x = 17M , m v = 55M , and K x = 283 km/s, still 
too high to be acceptable. Less than half of the stellar 
surface is sufficiently cool to give the absorption lines un- 
der study (e.g. due to sideway heating from scattered UV 
radiation in the strong accretion disk wind). This addi- 
tionally decreases the value of the actual radial velocity 
semi-amplitude. Fig. 1191 illustrates these considerations. 

So taking, for example, K v — 70 km/s and q = 0.3 
yields the optical star mass function f v = 0.46M Q , binary 
masses m x = 9Mq, m v = 30M Q and the optical star ra- 
dius R v ~ 4Oi?0. This radius is compatible with a typical 
bolometric luminosity of a 30M Q A5-A7 supergiant with 
T e ff ~ 8500 K. In this solution, K x ~ 233 km/s, larger 
than the measured value 175 km/s, but the true value of 
K x may be affected by the strong accretion disk wind. 

9. Discussion 

The interpretation of the binary parameters presented 
above is strongly based on the binary mass ratio q = 0.3 
as inferred from modeling the INTEGRAL observations of 
the hard X-ray eclipse with account of the observed ampli- 
tude (~ 4) of the hard X-ray precessional variability of the 
system. The model we use (an optical star filling its Roche 
lobe + thick accretion disk with a hot corona) is simplistic 
and cannot perfectly reproduce the apparent asymmetric 
shape of the X-ray eclipse. In the real situation, an ad- 
ditional X-ray absorption by outflowing gaseous streams, 
extended stellar envelope or wind seems quite plausible. If 
we interpret the complex shape of the hard X-ray egress 
as being due to additional variable absorption effects, we 
obtain the geometry of the hard X-ray emitting region in 
the form of an extended oblate corona over the accretion 
disk. However, a straightforward interpretation of the joint 
RXTE+INTEGRAL spectrum of the source observed in 
March 2004 is also possible by a single thermal emission 
with a temperature of ~ 30 keV. No additional emitting 
region is required to fit this broadband X-ray spectrum. 

There are two possible solutions. First, we could try 
to describe the observed hard X-ray spectrum by a more 
complex model involving two components, a hot corona 
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+ cooler thin jet radiating in the standard X-ray band. 
Second, we could try to accept the thin cooling jet inter- 
pretation of the X-ray spectrum, but then would have to 
admit a different shape of the X-ray eclipse. Our modeling 
shows that a thin short jet could reproduce the observed 
amplitude (~ 4) of hard X-ray precessional variability in 
SS433. In that case, the shape of the X-ray eclipse must 
be more sharp, with shorter ingress/egress times than ac- 
tually observed (cf. Fig. . If we look at the data, such 
a shape can indeed be found at the egress phase of the 
eclipse (see Fig. The (apparently longer) ingress time 
could be due to additional (apart from the opaque star 
body) absorption of hard X-ray emission by an optically 
thick gaseous stream outflowing from the optical star. 
Apparently the same smooth ingress and egress shape of 
the eclipse in softer X-rays (the Ginga and ASCA data) 
in that case could be due to the soft X-ray emitting region 
located further upstream of the jets. This interpretation 
also explains the shallower soft X-ray eclipse. 

Clearly, in this situation we need more observations 
of the hard X-ray eclipse to confirm the asymmetric 
ingress/egress eclipse shape and flatness of its bottom. 
We hope to do this in future INTEGRAL observations 
of SS433. 

10. Conclusions 

The results of the INTEGRAL observations of the peculiar 
binary system SS433 in 2003-2004 can be summarized as 
follows. 

1. Hard X-ray emission up to ~ 100 keV was dis- 
covered from this superaccreting microquasar with un- 
eclipsed hard X-ray luminosity close to the maximum 
opening precession phase L x (18 — 60keV) ~ 4 x 10 35 erg/s, 
L x (60- 120keV) - 2 x 10 35 erg/s (assuming the 5 kpc dis- 
tance to SS433), which is about 10% of the soft X-ray jet 
luminosity. 

2. Persistent precessional variability in hard X-rays was 
discovered with the maximum to minimum flux ratio ~ 4, 
which is twice as large as in the softer X-ray band. 

3. The observed hard X-ray orbital eclipse is found 
to be in phase with the optical and near infrared eclipses. 
The hard X-ray eclipse is observed to be at least two times 
deeper than the soft X-ray eclipse. The width of the X-ray 
eclipse increases with energy, which is opposite to what is 
observed in classical X-ray binary systems. 

4. The broadband X-ray spectrum 2-100 keV of SS433 
simultaneously obtained by INTEGRAL and RXTE in 
March 2004 when the source was in its flaring state can 
be fitted by bremsstrahlung emission from optically thin 
plasma with a temperature kT w 30 keV. 

5. Optical spectroscopic observations of SS433 on the 
SAO 6-m telescope which were performed in the frame- 
work of the multiwavelength INTEGRAL campaign of the 
source confirmed the spectral class of the optical star to 
be A5-A7I. The radial velocity curve of the optical com- 
ponent was obtained with a semi-amplitude K v = 132 ± 9 
km/s near the maximum disk opening precession phase 



(T3). The spectroscopic data revealed a heating effect of 
the hemisphere of the optical star illuminated by the ac- 
cretion disk. 

6. From the analysis of the hard X-ray eclipse and pre- 
cession variability in SS433 we estimate the mass ratio in 
this binary system to be q ~ 0.3. This mass ratio and 
radial velocity measurements corrected for the heating ef- 
fect enable us to evaluate the masses of the optical and 
compact star in SS433 to be M v w 3OM , M x w 9M , re- 
spectively. This finding lends further support to the pres- 
ence of a black hole in this peculiar superaccreting galactic 
microquasar. 
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